Black Hole Chromosphere at the LHC 
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If the scale of quantum gravity is near a TeV, black holes will be copiously produced at the LHC. 
In this work we study the main properties of the light descendants of these black holes. We show 
that the emitted partons are closely spaced outside the horizon, and hence they do not fragment into 
hadrons in vacuum but more likely into a kind of quark-gluon plasma. Consequently, the thermal 
emission occurs far from the horizon, at a temperature characteristic of the QCD scale. We analyze 
the energy spectrum of the particles emerging from the "chromosphere", and find that the hard 
hadronic jets are almost entirely suppressed. They are replaced by an isotropic distribution of soft 
photons and hadrons, with hundreds of particles in the GeV range. This provides a new distinctive 
signature for black hole events at LHC. 



General idea: In the mid-70 's Hawking 1] pointed out 
that a black hole (BH) emits thermal radiation as if it 
were a black body, with a temperature Tbh inversely pro- 
portional to its mass Mbh- More recently. Heckler 2] 
noted that for Tbh > Aqcd, the emitted partons do not 
immediatly fragment into hadrons, but propagate away 
from the BH through a dense quark-gluon plasma los- 
ing energy via QCD bremsstrahlung and pair produc- 
tion, and forming a nearly thermal chromosphere. As- 
suming that TeV scale gravity is realized in nature J|J , in 
this work we first show that such chromosphere would be 
shielding the BHs to be produced at the Large Hadron 
Collider (LHC) 0,13, and then we estimate the emer- 
gent thermal photon spectrum from tt" decay. The en- 
suing discussion will be framed in the usual context of 
TeV-scale gravity scenarios, in which the standard model 
(SM) fields are confined to the brane and only gravity 
spills into the n internal dimensions of the universe. We 
also rely on the probe brane approximation, viz., the only 
effect of the brane field is to bind the BH to the brane. 
This will be an adequate approximation provided Mbh 
is well above the brane tension, which is presumably of 
the order of but smaller than the fundamental Planck 
scale, Mb. Moreover, we assume the BH can be treated 
as a flat A-\-n dimensional object. This assumption is 
valid for both extra dimensions that are larger than the 
Schwarzschild radius |a, 
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and for warped extra dimensions where the horizon ra- 
dius is small compared to the curvature scale of the ge- 
ometry associated with the warped subspace . 

Hawking evaporation: According to the Thome's hoop 
conjecture 0, an event horizon forms when and only 
when a mass AfsH is compacted into a region whose 
circunference in every direction is less than 27rrs. The 
strong gravitational fields around the BH induce sponta- 
neous creation of pairs near the event horizon |lj . While 
the particle with positive energy can escape to infinity. 



the one with negative energy has to tunnel through the 
horizon into the BH where there are particle states with 
negative energy with respect to infinity p. As the BHs 
radiate, they lose mass and so will eventually evaporate 
completely and disappear. The evaporation is generally 
regarded as thermal in character, with a temperature 
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and an entropy S = 4 tt Mbh ''^/(''^ + 2). However, the 
BH produces an effective potential barrier in the neigh- 
borhood of the horizon that backscatters part of the out- 
going radiation, modifing the blackbody spectrum. The 
BH absorption cross section, CTs (a.k.a. the greybody 
factor), depends upon the spin of the emitted particles 
s, their energy Q, and the mass of the BH Mbh 
At high frequencies {Qrs ^ 1) the greybody factor for 
each kind of particle must approach the geometrical op- 
tics limit. For illustrative simplicity, we adopt through- 
out this work the geometric optics approximation, and 
following ^2 1 we conveniently write the greybody factor 
as a dimensionless constant, Tg = ag/A^, normalized to 
the BH surface area 
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seen by the SM fields (Fs^o = 1, ~ 2/3, and 

Fs=i ~ 1/4 01). The prevalent energies of the decay 
quanta are ~ Tbh ~ l/'^s: resulting in s-wave dominance 
of the final state. This implies that the BH decays with 
equal probability to a particle on the brane and in the ex- 
tra dimensions jl3l |. Therefore, the evaporation process 
is dominated by the large number of SM brane modes. 
Since the s-wave greybody factors for fermions |l j | dif- 
fer from the geometric ones by only about 10-20%, we 
continue to use the latter. 

A A+n BH then emits particles with initial total energy 
between [Q, Q -|~ dQ) at a rate 
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per degree of particle freedom i. Integration of Eq. Q 
leads to 
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where T{x) (C(a^)) is the Gamma (Riemann zeta) function 
and / = 1 (/ = 3/4) for bosons (fermions). Therefore, 
the BH emission rate is found to be 
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FIG. 1: Feynman diagram for the dominant contribution to 
gluon bremsstrahlung in a qq or qq collision. 
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for particles with s — 0, 1/2, 1, respectively. Note 
that the BHs to be produced at the LHC are per- 
fectly well defined resonances with mean lifetimes, r^^ ^ 



10 



-27 



s, and Tbh 



200 GeV. 

Since thermal fluctuations due to particle emission are 
small when the entropy 5 3> 1 fll and statistical fluc- 
tuations in the microcanonical ensemble are small for 
y/S 3> 1 3, we expect the above formulae to be an ad- 
equate approximation for Mbh ^ Mij. This condition 
inevitably breaks down during the last stages of the de- 
cay process. However, it is noteworthy that for BHs with 
initial masses well above Md most of the evaporation 
takes place within the semi-classical regime. In what fol- 
lows we require Mbh/Md ^ 5, corresponding to S" > 25. 
To obtain some quantitative estimates on the particle 
spectra, we hereafter set n = 6 and Mu = 1.3 TeV, a 
value consistent with current limits on the size of extra- 
dimensions 

Thermalization: The evaporation of the BH creates a 
radiation shell of radius and thickness t^^ , which prop- 
agates outward at the speed of light. Using the preceding 
equations we find that, already for Mbh/ Mo = 5, there 
are about 10 quarks and antiquarks in the shell. Before 
proceeding further, we examine under what conditions 
these particles can undergo sufficient dissipative interac- 
tions (bremsstrahlung and pair production) to initiate 
thermalization. 

At any time t after the end of evaporation, the shell 
has radius r = ct and thickness t^^ . A quark or antiquark 
in the shell will interact at a rate given by 
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where n is the qq effective density at time t, and cr^^^^ 
is the cross section for gluon bremsstrahlung in a qq or 
qq collision. The dominant diagram of this process at 



low momentum transfer is shown in Fig_ 
averaged total cross section is given by 
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where as is the QCD coupling constant, and A is an in- 
frared cutoff related to the off-shell momentum of the ex- 
changed gluon. Note that the exchanged (virtual) gluon 
can travel at most a distance ~ r during the interaction, 
so that the uncertainty principle provides a lower bound 
~ on the momentum transfer. Thus, at any time 
t = r/c, the cross section has a maximum value corre- 
sponding to A « r~^. Of course, A is absolutely bounded 
from below by Aqcd- 

In light of these considerations, we take the effective 
density to be fully volumetric 



so that a quark or antiquark interacts at a rate 
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For 
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we can choose A 



In order to 



remain perturbative, we find the cumulative number of 
interactions/quark for A ^ down to some value 3> 
AqcD- Integration of ((T^ with Nqq w 10 leads to 
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The fiducial value of as has been chosen appropriate to 
a momentum transfer scale of 9 GeV |l6j . correspond- 
ing to upsilon decay. Even for such a large value of 
A and g ~ 2 Tbh ^ 400 GeV, Eq. ^ yields « 3 in- 
teractions per quark; for A — 1.78 GeV, corresponding 
to T-decay [Hi, ~ 0.35, and Eq. ^ yields « 30 
interactions per quark. The creation of a gluon com- 
ponent in the plasma will lead to a significant increase 
in the bremsstrahlung cascade, because the gluon-gluon 
cross section is considerably larger than the one of quark- 
(anti)quark. We take this result as supportive of thermal- 
ization, and we proceed on that assumption. 
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Before interactions have led to a rapid increase in par- 
ticle numbers, the average separation at distance r of 
partons emitted by the BH is defined as 
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where Ci is the number of internal degrees of freedom of 
particle species i, i.e., 72 (16) for quarks and antiquarks 
(gluons). Evaluating this expression for our parameters, 
we find that (d)parton ~ 7 x 10"^GeV"^ < ^qcd- 
interactions cascade, the volume occupied by the partons 
expands as well. The temperature of the chromosphere 
decreases as the radius increases, and eventually reaches 
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at which point quarks and gluons hadronize. 

Soft thermal spectrum: The density of the outward- 
going plasma changes by a large factor within a mean 
free path, and so the particles never have enough time to 
fully thermalize. In order to describe the outer edge of 
the chromosphere we adopt here the heuristic treatment 
given by Heckler 0. The quark and gluon spectrum in 
the observer rest frame is obtained by boosting a ther- 
mal spectrum at = Aqcd with the Lorentz factor, 
7^j^ « (Tbh/Aqcd)^^'^, of the outer surface of the ther- 
mal chromosphere, i.e.. 



m 

dQ 



27r2 



dn{l- [3 COS 61) cos( 



exp 



X,Q(l-/3cos0) 



(-1) 



2s 



,(16) 



where the integration over the surface of the chromo- 
sphere of radius r^^ = 7^,, /Aqcd has been carried out. 
It is noteworthy that detailed analyses ^3 based on the 
transport equation yield results in very good agreement 
with this spectrum. 

The observed photon spectrum of the chromosphere 
is a convolution of the quark-gluon spectrum given in 
Eq. lUni) with the pion fragmentation function and the 
Lorentz-boosted spectrum from tt*^ decay; namely. 
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where Eq = E^ + rn^/AEj. The number of photons with 
energy Ery produced by a pion propagating with velocity 
P and decaying isotropically in its rest frame is 

dg^^{E^) _ 2 _ 2 



~dE^ 



im^P {El-miyn' ^^^^ 

with 7 = (1 — /3^)~^/^. Finally, the pion spectrum can 
be expressed in the form .18.] 
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FIG. 2: Photon emission spectra of a 10-dimensional 
Schwarzschild BH, assuming Md ~ 1.3 TeV and 
Mbh/Md = 5. This implies, Tbh ^ 200 GeV and 
Tgjj ~ 4 X 10~^^ s. The solid line denotes the spectrum from 
quark fragmentation and subsequent tt" decay, whereas the 
dashed line stands for direct photon emission from the BH. 



where the sum is over relevant species in the plasma. The 
precise nature of the fragmentation process is unknown. 
As in our previous analysis of BH decay '^^ adopt 
the quark hadron fragmentation function originally 
suggested by Hill ^2Qj 
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that is consistent with the so-called "leading-log QCD" 
behavior and seems to reproduce quite well the multiplic- 
ity growth as seen in collider experiments {z = Et^/Q). 

Using Eq. H17|) one can easily estimate the mean life- 
time of the chromosphere 
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where AEj is the full width at half maximum (max) 
of the dN/dE^ distribution. The spectrum peaks at an 



energy of E^ 



100 MeV, and thus for our choice of 



parameters one obtains Tch ~ 2.2 x 10^ . 

In Fig. 13 we show the emergent photon spectrum from 
the chromosphere hadronization and subsequent 7r° de- 
cay, together with that of direct photon emission from the 
BH. The photon emission rate from the chromosphere is 
found to be 0] 
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Therefore, for a BH with w 4 x 10"^'^ s and Tbh ~ 
200 GeV, one expects ~ 7 x 10"* photons as end prod- 
ucts of the quark-gluon plasma. The nominal Compact 
Muon Solenoid (CMS) design anticipates the dynamic 
range for the electromagnetic calorimeter readout will ex- 
tend from the noise floor of 25 MeV up to 2 or 2.5 TeV, 
whereas the hadron calorimeter readout will range from 
20 MeV to 2 TeV j2l|. Consequently, detection of chro- 
mosphere's radiation could become observationally fea- 
sible. Moreover, a tt^ spectrum with roughly the same 
shape, but shifted in energy by a factor of 2, would trig- 
ger the hadron calorimeter. 

In general, one expects a statistically significant si gna l 
at the LHC from both pp j3| and PbPb collisions |22l |. 
In particular, for the above parameters and assuming a 
geometric cross section for the colliding partons, about 
10^ BH events with zero net transverse momentum px 
would be produced with an integrated luminosity of 
100 fb~^ Proposed modifications, such as BH pro- 
duced with high pt [l^l and recoil effects from graviton 
emission into the bulk [3| do not qualitatively change 
this estimate. Criticisms of the absorptive black disc 
scattering amplitude, which center on the exponential su- 
pression of transitions involving a (few-particle) quantum 
state to a (many-particle) semiclassical state [2^, have 
been addressed in [2^. However, it is noteworthy that 
even if one includes this supression, the BH production 
rate will still be quite large at the LHC for over almost 
all of interesting parameter space Specifically, for 

our fiducial values about 10'' BHs (either at rest or trav- 
elling along the beam pipe) would be produced with an 
integrated luminosity of 100 fb""'. 

Note that there is a mild dependence of Tbh with 
the initial BH mass. Consequently, most of the BHs 
{Mbu/Mu > 5) to be produced at the LHC have temper- 
atures « 200 GeV. The spectrum from quark fragmen- 
tation and subsequent tt^ decay given in Fig. El is thus 
quite general. 

Summary and final remarks: The results of this work 
indicate that the signals for BH production delineated 
by Giddings-Thomas and Dimopoulos-Landsberg 
(large transverse momentum, hard electromagnetic com- 
ponent) can be strikingly augmented with the observa- 
tion of the soft thermal isotropic hadronic and gamma 
ray component, which includes significant photon and 
charged pion counting rates in the vicinity of ~ 1 GeV. As 
can be seen from Fig. one expects about 100 photons 
(as well as about 200 charged pions) per black hole event 
with energies in the GeV range. A hard hadronic compo- 
nent is largely absent. This is much more dramatic than 



BH 



the absence of hadronic jets with pT > r^ 
According to calorimeter design guidelines for the CMS 
detector, the observation of this component should be 
feasible. 

TeV-scale BHs can also be produced in cosmic ray col- 
lisions li^, and (of course) one expects the quarks and 



gluons to be at first freely streaming away from the hori- 
zon. However, contrary to collider experiments, the BH 
is produced with large momentum in the lab system, and 
its decay products are swept forward with large, O(IO^), 
Lorentz factors. This minimizes the effect of the chro- 
mosphere on the spectra of the emergent particles ana- 
lyzed in Ref. In particular, a 100 MeV pion boosted 
by 10^ will still have a mean decay length to muons of 
~ 8000 km, much larger than 1 km, which is the esti- 
mated mean interaction length of a pion of 10^ GeV [3nj |. 
Thus, an extensive air shower will still occur. The final 
detectability of the shower will differ little from the situ- 
ation when the decay particles have energies ^ 100 GeV. 
If anything, the chromosphere will accelerate the muon 
production, for reasons similar to the enhancement of 
the muon component in a shower initiated by a complex 
nucleus. 

Our analysis impacts directly on the identification of 
the top quark among the BH subproducts j^^l • The decay 
rate of the t-quark in the SM is approximately 

~ t|- "^t (Gf/V2) ~ 0.5 GeV , (23) 

yielding a lifetime T( ^ 2 GeV^^. Here, (m^i/) denotes 
the mass of the top (W) and Gp is the Fermi coupling 

constant. Now, since rch ~ ^bh /^qcd ~ 150 GeV~^ 
(and the effects of time dilation are negligible), top 
quarks emitted as part of the Hawking radiation cannot 
escape the chromosphere even in the absence of collisions. 
Diffusive effects will increase the time spent in the chro- 
mosphere, and thus strengthen this result. Therefore, 
top quark identification depends on the leptonic activity 
resulting from their decay. We also note that both the 
W and Z decay inside the chromosphere. On the other 
hand, the lifetime for a Higgs boson decaying into bb is 
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where mt {rriH) is the mass of the 6-quark (Higgs) and 
V = 246 GeV. In two-Higgs models (such as supersym- 
metry) , tan (3 is the ratio of the vevs and a is a mixing 
angle in the neutral Higgs sector. For the SM one-Higgs 
case, the factor cos (3/ sin a is set equal to 1. Thus, in the 
SM, the Higgs seems likely to escape the chromosphere 
and generate a hard ietty component in small percentage 
of the BH decays [S^. In the two-Higgs model, escape 
from the chromosphere will depend on the mixing angles. 
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